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SUMMARY 


Charts are presented to provide a convenient means for obtaining 
the derivatives of rotor resultant force, rotor pitching moment about 
the helicopter center of gravity, and rotor torque with respect to rotor 
angle of attack, forward speed, rotor speed, and collective pitch. The 
analysis is made for untwisted, untapered blades which have flapping 
hinges located at the rotor shaft. The charts can, however, be used 
without serious error for blades with normal values of taper and twist 
and means are indicated for obtaining the additional terms necessary when 
the charts are applied to rotors having moderate flapping-hinge offsets. 
These charts do not rely upon the assumption (often used for convenience) 
that the rotor resultant-force vector is perpendicular to the tip- path 
plane. This assumption is shown to lead often to grossly incorrect 
longitudinal- stability derivatives . 

The method of using the charts to obtain the various derivatives is 
explained and illustrated. The use of the charts in the study of the 
general stability characteristics of the helicopter rotor is also 
illustrated. 


INTRODUCTION 


During the past few years, the factors that affect the performance 
of the helicopter have become fairly well understood and, as a result, 
the performance of the helicopter has been improved appreciably. Refer- 
ences 1 and 2 indicate, however, that the flying qualities of helicopters 
at the present time are still, in general, inferior to those of airplanes. 
Reference 1 states that future Navy specifications will incorporate 
stringent flying-qualities requirements for helicopters. According to 
reference 2, the longitudinal flying qualities in forward flight are most 
in need of improvement. Reference 3 presents tentative longitudinal 
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f. lying-qualities requirements based on the normal-acceleration charac- 
teristics during a pull-up maneuver in forward flight. This reference 
indicates that a helicopter that meets these requirements will be much 
safer and less fatiguing to the pilot than one which does not. 

In order to determine whether a prospective helicopter design will 
meet the longitudinal flying-qualities requirements, stability calcu- 
lations for a pull-up maneuver must be made. Such a calculation involves 
setting up and solving simultaneously the differential equations of 
motion of the helicopter. In order to set up these equations, the sta- 
bility derivatives of the helicopter and, therefore, of the rotor must 
be known. The stability derivatives are also useful for prediction of 
other stability and control characteristics and for the study of stabi- 
lizing devices and autopilots. In addition, these derivatives provide a 
convenient basis for estimating load factors for arbitrary gust values. 

The rotor theory of references 1|. and £ has been used for performance 
and blade-motion calculations. These calculations have compared favor- 
ably with flight test results (see, for example, references 6 and 7); 
therefore, experimental checks are expected to show the same theory suit- 
able for stability calculations. The results of reference 3 show that 
this theory can be used for analytical studies of the ability of the 
helicopter to meet the tentative longitudinal flying-qualities require- 
ments of that reference. This theory, however, is presented in a form 
convenient for performance and blade-motion studies but is not particularly 
convenient for stability studies. Thus, a rework of the theory in order 
to present it in a form more convenient for stability studies was con- 
sidered desirable. For example, in estimating horizontal-tail area needed 
to offset the angle-of-attack instability of a helicopter rotor, it was 
found that, although the changes in magnitude of the rotor resultant-force 
vector with angle of attack could readily be calculated, estimation of 
the changes in vector longitudinal tilt with angle of attack to the 
necessary degree of accuracy was relatively difficult. Approximation of 
the vector longitudinal tilt by the easily calculated longitudinal flap- 
ping a^ appeared excessively inaccurate for the purpose. The basic 
theory was accordingly used to prepare the charts presented herein from 
which the values of the parameter a', which represents vector longitudinal 
tilt, could be read directly. For convenience, values of torque coef- 
ficient are also obtainable from these charts, and additional charts from 
which variations in magnitude of the resultant-force vector can be read 
directly are also included. 

The charts presented herein can be conveniently used to obtain most 
of the rotor longitudinal-stability derivatives. Those derivatives which 
depend upon pitching velocity are not included in the present paper; 
however, the damping moment due to a steady pitching velocity, which is 
the most significant pitching-velocity derivative, can be obtained from 
reference 8. 



NACA TN 2309 


3 


SYMBOLS 

Physical Quantities 


b 

r 

R 

c 


CT 

9 


*1 

P 

y 

6 3 


number of blades per rotor 
radial distance to blade element, feet 
blade radius, feet 
blade-section chord, feet 


equivalent blade chord (on thrust basis), 



rotor solidity 



blade-section pitch angle; angle between line of zero lift of 
blade section and plane perpendicular to axis of no feathering, 
radians 


mass moment of inertia of blade about flapping hinge, slug-feet 2 


mass density of air, slugs per cubic foot 


mass constant of rotor blade, 


to mass forces 


(cpaR^/l]_) 


expresses ratio of air forces 


angle in plane of rotation between perpendicular to blade-span 
axis and flapping-hinge axis, positive when an increase in 
flapping produces a decrease in blade pitch 


Air-Flow Parameters 

V true airspeed of helicopter along flight path, feet per second 

V h horizontal component of true airspeed of helicopter, feet 

per second 


vertical component of true airspeed of helicopter, feet 
per second 

rotor angular velocity, radians per second 

rotor angle of attack; angle between flight path and plane 
perpendicular to axis of no feathering, positive when axis 
is pointing rearward, radians 

tip-speed ratio 

induced inflow velocity at rotor, always positive, feet 
per second 


V cos a 
^ £2R 


assumed equal to — 

H m 


inflow ratio f V . 5i “ - Q ~ - X 
\ fiR 

blade azimuth angle measured from downwind position in direc- 
tion of rotation, radians 

component at blade element of resultant velocity perpendicular 
to blade-span axis and to axis of no feathering, feet 
per second 


component at blade element of resultant velocity perpendicular 
both to blade-span axis and Uip, feet per second 

inflow angle at blade element in plane perpendicular to blade- 
span axis, radians ^tan -1 

blade-element angle of attack, measured from line of zero lift, 
radians (9 + 0) 



Aerodynamic Characteristics 
section profile drag coefficient 

slope of curve of section lift coefficient against section 
angle of attack, per radian 
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L rotor lift, pounds 

T rotor thrust, component of rotor resultant force parallel to 

axis of no feathering, pounds 

Q rotor- shaft torque, pound-feet 


rotor thrust coefficient 


,nR 2 p(fiR) 2 / 


Cq rotor-shaft torque coefficient 


> nR 2 p(nR) 2 Ry 



Performance Parameters 
rotor profile drag-lift ratio 



rotor induced drag-lift ratio 



parasite drag of helicopter components other than lifting 
rotors divided by rotor lift 



drag-lift ratio representing angle of climb or glide, positive 
in climb ( V v /V h ) 




P/L 


rotor drag-lift ratio; ratio of equivalent drag of rotor to 
rotor lift ((§) o * (2).j 

component of rotor resultant force along flight path (that is, 
useful component of rotor resultant force) divided by rotor 

lm ((e)p + (£)o) 

shaft- power parameter, where P is equal to rotor-shaft power 
divided by velocity along flight path and is therefore also 
equal to drag force that could be overcome by shaft power 

at flight velocity \ 
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angle of climb 




Rotor-Blade Motion 

(3 blade flapping angle at particular azimuth position, radians 

a q constant term in Fourier series that expresses (3; therefore, 

rotor coning angle 

a-^ coefficient of -cos in expression for J3; therefore, 

longitudinal tilt of rotor cone 

b-^ coefficient of -sin in expression for (3; therefore, lateral 

tilt of rotor cone 


Stability 

a' projection of angle between rotor force vector and axis of no 

feathering in plane containing flight path and axis of no 
feathering 

q helicopter pitching velocity, radians per second 

A increment 

t a b functions of p given in reference Uj subscripts a and b 

" ’’ represent numbers used to identify a particular function 

C parameter defined by equation (A10) in appendix A 

STABILITY DERIVATIVES NEEDED FOR PULL-UP ANALYSIS 


The stability derivatives needed for a pull-up analysis are now 
discussed, not only because of the significance of such an analysis but 
also because most of the derivatives needed for other longitudinal- 
stability studies are included in this relatively complex case. The 
actual use of these derivatives in a pull-up analysis is not discussed 
herein. 

Rotor forces and blade flapping are affected by five independent 
variables: rotor angle of attack a, forward speed V, rotational 

speed fl, collective pitch 9, and pitching velocity q. The helicopter 
is affected by both the magnitude and direction of the rotor resultant 
force and by the magnitude of the rotor torque. (The rotor resultant-force 
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vector is assumed, to pass through the rotor hub inasmuch as present-day- 
conventional rotors are designed with near-zero airfoil pitching moments 
and with coinciding aerodynamic centers and chordwise centers of gravity. ) 
If the flapping hinges are offset from the rotor shaft or cocked, the 
helicopter will also be affected by longitudinal and lateral flapping and 
coning-angle changes. The magnitude and longitudinal direction of the 
rotor resultant force depend upon Cj/a and a', respectively, and the 

magnitude of the rotor torque depends upon Cq/ct. Calculation of the 

helicopter stability derivatives needed for a pull-up analysis therefore 
requires knowledge of the following rotor derivatives: 

8(0 t /o) 3(C T /a) d(C T /a) 9 ( c t/ ct ) d(c T /a 

3a 3V 3fi 30 3q 

3a* 3a 1 3a 1 3a 1 3a 1 

3a 3V 3ft 39 3q 

8(Cq/o) 3(Cq/o) 3(Gq/o) 9(Cq/o) S(Cq/o) 

3a 3V 3 0 30 3q 


3a]_ 

3a]_ 

3a^ 

3a^ 

3a-^ 

3a 

3V~ 

3ft~ 

30“ 

3q 

3a 0 

3a 0 

3a 0 

3a 0 

3a 0 

3a 

3V 

3ft 

30 

3q 

3b x 

3b x 

db 1 

3b 

3b-L 

3a 

3V~ 

dn 

30~ 

3q 


By definition, each partial derivative is obtained by assuming the other 
four independent variables to be constant. 

All the preceding Cip/o, a', and Cq/c derivatives, except those 

with respect to q, are obtainable from the charts in this paper. The 
aQ, a-^, and b-^ derivatives, except those with respect to q, can be 

obtained from the equations of reference U. Although reference U is 
restricted to the case of flapping hinges on the rotor shaft, the equa- 
tions are believed to be reasonably accurate for small values of flapping- 
hinge offset. The derivatives with respect to q are not included in 
the present paper; however, the changes in a' and a^_ due to a steady 

pitching velocity, which are the most significant pitching-velocity 
derivatives, can be obtained from reference 8. 
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For specific helicopter configurations, many of these derivatives 
may be expected to be insignificant, but no specific investigation of 
this possibility has yet been made. 

For certain longitudinal-stability calculations, other rotor aero- 
dynamic derivatives are desired. For example, in order to determine the 
rotor contribution to the variation of stick position with speed at 
fixed pitch, thrust, and rotational speed, the derivative da'/cSV at 
constant Cp/o (instead of at constant a) is needed. This a' deriva- 
tive can also be obtained from the charts in this paper. Similarly, 
da' /dQ and da'/dQ, as well as the derivatives of Cq/ct with respect 

to V, fl, and 0, can be obtained at constant Cp/d from the charts 
presented herein. 


ASSUMPTIONS OF THEORY 


In order to limit the labor involved in the theoretical derivations 
necessary to prepare the charts of this paper, various assumptions were 
made. For the readers' convenience, the more significant assumptions 
carried over from references U and 9 as well as those made herein are 
discussed. 


Rotor Physical Characteristics 

The rotor is assumed to have freely flapping rectangular untwisted 
blades with mass factor y of 15 and the flapping hinge is assumed to 
be on the rotor shaft perpendicular to the blade-span axis. The charts, 
however, are considered applicable to rotors with blades having a range 
of y from 0 to 25 , moderate amounts of twist or taper, and moderate 
amounts of flapping-hinge offset or angularity ^6^ angle 

The consideration that the charts are applicable to rotors having 
values of y ranging from 0 to 25 is based on the fact that the theory 
is shown in reference U to be applicable to this range of y when a 
value of y = 15 is used. The fact that the theory is satisfactory over 
this wide range of y, and therefore for a variety of values of coning 
angle, lateral flapping angle, and of higher harmonics of flapping, 
suggested that it might be applicable to see-saw rotor systems. Further 
investigation of the problem indicated that the charts presented herein 
should, in general, give sufficiently accurate results for see-saw rotor 
systems . 
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The charts are considered applicable to rotors having moderate 
amounts of taper because the satisfactory comparison between calculated 
and measured rotor performance and blade motion previously referred to 
was obtained with blades which actually had a moderate amount of taper. 

The taper was accounted for theoretically by basing the solidity a and 
the blade mass factor y on an equivalent blade chord. (See "Symbols.'') 

The consideration that the charts are also applicable to rotors 
having moderate amounts of twist is based on a comparison made between 
the charts presented herein and calculations for a rotor with -8° of 
linear twist. This comparison showed that the derivatives are not 

significantly affected. The absolute values of a 1 given by the charts 
are nowhere in error by more than 1° and rotor instability with angle of 
attack is indicated to be somewhat high, which is conservative. The 
other a' derivatives and the Cq/o derivatives for the rotor with 

-8° twist are given by the charts to a satisfactory degree of engineering 
accuracy. 

The applicability of the charts to rotors with moderate amounts of 
flapping-hinge offset is logically justified as follows. Small amounts 
of flapping-hinge offset (less than 5 percent R) are not considered 
likely to affect rotor performance. Thus, inasmuch as the charts presented 
herein are derived from performance equations, they should be applicable 
without excessive error to rotors having flapping-hinge offsets of less 
than 3 percent of the rotor radius. For such rotors, however, the effect 
of longitudinal flapping with respect to the shaft in producing additional 
pitching moments due to the centrifugal force in the blades must be taken 
into account. 

The charts can be applied to rotors having angular flapping hinges 
(&3 angle) if the resulting effect of lateral flapping on the longitudinal 
position of the axis of no feathering is taken into account. (The axis of 
no feathering is the axis about which no first harmonic pitch variation 
occurs. See appendix of reference 7 for discussion of its use and 
significance . ) 


Section Aerodynamic Characteristics 

The equation representing the blade-section profile drag coefficient 

is 

c do = 0.0087 - 0 . 02l6a r + 0.1|00a r 2 

Comparisons between this theoretical profile drag polar and statically 
measured experimental polars for several typical airfoils are given in 
figure 1 of reference 5. This theoretical polar can also be compared with 
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statically measured drag polar data presented in reference 9 . The sec- 
tion is assumed to have a lift-curve slope of 5.73 per radian. This 
value can be compared with statically measured lift curves of various 
rotor-blade airfoil sections which are also presented in reference 9* 

These lift and drag section characteristics are assumed to be con- 
stant over the rotor disk; that is, effects of radial velocities and 
variations around the disk in Reynolds number and Mach number are 
neglected. 


Rotor Aerodynamic Characteristics 

In the calculation of rotor aerodynamic characteristics from the 
blade-section characteristics, the inflow velocity is assumed to be 

uniform and 9, 0, a , (3, a, tan“-*-(2.) } and a* are assumed to be 

\Vu 

small angles so that the sine is equal to the angle and the cosine is 
unity. The charts are limited to values of p. between 0.15 and 0.50 
where these assumptions are believed to be valid. The most questionable 
assumption is that of uniform inflow velocity. This assumption is con- 
sidered reasonable for calculating longitudinal- stability derivatives in 
view of the previously referred to satisfactory comparisons between rotor 
theory based on this assumption and both measured longitudinal flapping 
and rotor performance. 

The loss in lift at the tip of the blade is approximated in cal- 
culating rotor aerodynamic characteristics by assuming that the outer 
3 percent of the blade radius is ineffective in producing lift but does 
produce drag. 

In calculations of rotor thrust, the contribution of blade-section 
drag is omitted as being negligible. 

LIMITS OF VALIDITY OF THEORY 
Stalling Limits 

As was explained in references U and 5, the accuracy of the theory 
becomes doubtful when a stalling angle of attack is encountered at ^ = 270° 
by either a section at the blade tip or by an inboard section with relative 
velocity equal to four-tenths the rotational tip speed. The reasons for 
this doubtful accuracy are that the theory fails to take account of the 
nonlinearity of the lift curve and the increased drag rise at the stall. 
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The charts of this paper, therefore, follow the procedure of the 
charts of reference 5 in that locus lines for 12° and 16° angles of 
attack are included. The locus lines for the conditions for which a 
blade element at an azimuth angle of 270° with a relative velocity upRR 
equal to four-tenths the rotational tip speed reaches a specified angle 
of attack are designated by the symbol a ^ UT=0 ^)(270°)* Sirailarl y> the 


locus lines for the conditions for which the blade tip at an azimuth 
angle of 270° reaches a specified angle of attack are designated by the 
symbol o)( 270 °) - Inasmuch as a rotor with twisted blades has approxi- 

mately the same pitch at the three-quarter radius as a rotor with untwisted 
blades for the same flight condition, the specified angles on the tip locus 
lines should be reduced by approximately one-fourth of the blade twist 
value when using the charts for twisted blades. A similar correction that 
depends upon p can be calculated for the specified angles on the locus 
lines for Up = 0.1;. 


As pointed out in several previous MCA papers (see, for example, 
reference 5) , the optimum flight condition from considerations of heli- 
copter performance is approximately that at which stalling just begins. 

A helicopter designed to fly near this optimum performance condition 
therefore always enters the stalled region during an appreciable pull-up 
or upward gust. It would therefore be desirable to determine the 
necessary empirical correction, if any, to the stalled part of the charts 
presented herein. This correction, however, is believed to be beyond 
the scope of the present paper. 


Mach Number Limits 

Although the section lift and drag characteristics vary with Mach 
number, on the basis of performance prediction experience, it is believed 
that such changes would not significantly affect the validity of the 
charts presented herein (except perhaps fig. 1(a)) at least until the. 

Mach number for lift and drag divergence is reached. It is possible that 
these charts may be valid somewhat beyond this point because of a tip- 
relief effect and other factors, but further discussion of this point is 
believed to be beyond the scope of this paper. 

The quantities given in figure 1(a) are more directly dependent on 
lift-curve slope and, for helicopters with unusually high tip speeds, the 
effect of Mach number may be significant even before the lift divergence 
is reached. For such designs, the average section lift-curve slope may 
be appreciably higher than the value of 5.73 used herein. If such is 
the case, equations given in appendix A may be used in place of figure 1(a) 
for greater accuracy. 


Rotor Angle-of -Attack Limits 

Inasmuch as the term cos a occurs in the accurate expression 
for p, a change in rotor angle of attack while Constant forward speed 



12 


NACA TN 2309 


and rotor speed are maintained causes a small change in p. This change 
in p in turn causes a small correction to the derivatives with respect 
to a as calculated with cos a assumed to equal 1-0. If desired, this 
correction is obtained as follows: 

V cos a 

^ dr 

Differentiating this expression yields 


Therefore, 



const. 


V,Q 


V sin a 
flR 


= -p tan a 



const. 


V 



const. 



tan a — 
dp 


( 1 ) 


Calculations indicate that the rotor angle of attack is usually small 
enough to make this correction term insignificant. 


Need for Experimental Verification 

Because of the various assumptions used in the theory on which the 
charts presented herein are based, adequate experimental verification is 
necessary before these charts can be used with complete confidence. Only 
approximate and incomplete experimental checks are as yet available. 
Nevertheless, because of the satisfactory correlation of experiment and 
theory for rotor performance and blade motion previously referred to and 
because of the satisfactory comparison between measured and predicted 
pull-up time histories shown in reference 3, it is believed that these 
charts can be used with some confidence even before any further experi- 
mental verification is obtained. 


DERIVATION OF CHARTS 
Crj, jo Derivatives 


The Cijt jo 
of reference b. 


derivatives are obtained by combining equations (6) and (7) 
which (with the blade-twist term dropped) are, respectively, 


2C t 

— = * 3 , 1 * + * 3,29 
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and 


tan a = ^ + 

* 2n(x 2 * ^ 2 ) 

(The symbols such as t^ ]_ refer to tabulated constants in reference U. 

The subscript designates the table number and line. For a rotor with 
twisted blades, 9 is approximately the pitch at the three-quarter radius.) 



The manner in which these equations are combined and processed is 
described in appendix A. As explained in that appendix, typical values 
of X and are employed in obtaining figures 1(a) and 1(b) (X = -0.05, 

C T = 0.005). Calculations indicate that, for p ^ 0.15, this approximation 

should seldom involve any significant error but, for special cases, the 
equations given in appendix A can be employed. 

Derivatives of C^/a with respect to 9 and a.- The deriva- 

d(C T /a) 3 (c T /a) 

tive — j — - — for constant p and a and the derivative for 

39 3a 

constant p and 9 are plotted against p in figure 1(a) for values 
of p from 0.15 to 0.50 and for a = 0.03, 0.06, and 0.09. 


Derivative of (Cy/o'j 
for constant a and 9 


with respect to p.- The derivative 
is given by the following equation: 


9 (V°) 

3p 


S(C T / a ) 

3p 


k-j9 + k2 — + kjC<p9^ + k^ 


Cm^9 Crp^ 

— — + k^ -~ 

a 5 a 2 


(2) 


where the k quantities are plotted against p in figure 1(b) for 
values of p from 0.15 to 0.50 and for a = 0.03, 0.06, and 0.09. 


a 1 


and 



Derivatives 


Figure 2 shows that the longitudinal angle between the rotor 
resultant-force vector and the axis of no feathering is 




# 


arc tan ( : 

r) + a 
w u J 


( 3 ) 



Ill 
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Therefore , use of the trigonometric formula for the tangent of a 
sum of two angles gives 


tan a 1 


(r 


+ tan a 


u 


1 


-(? 


(U) 


tan a 


u 


Substituting equation (7) of reference h and equation (l) of 
reference 5 gives 


tan a* 


(R) + (D\ _P_X_ 

W 0 \LJ i L tx 


+ X 2 


1 - ( ^ tan a 

■ L Ai 


( 5 ) 


According to the second term of equation (ll-U) of reference 10, 


Thus, 


D\ = C T 
1 2(J.\1 + X 2 



(6) 


Substituting for X from equation (6) of reference U (dropping the blade- 
twist term) results in the following equation: 



2Crn 

to 2 0 

a a 


(it 


hL 



tan a 1 


tan a 


(7) 
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Assuming that a', a, and arc tan ^ are small angles, a' can 
be expressed in degrees as follows: 


i'° = 57.3 


In - i 


t 3 2 e ° - 57.3 
+ -hi 

t 3,l^ 


2C, 


T 

aa 


( 8 ) 


The angle a'° can be calculated from equation (8) as a function 
of [i , C t /ct, and 9 by first calculating (D/L) 0 and P/L as 

functions of p, Cr^/o, and 9 from the equations of reference l; as was 

done to obtain the charts of reference 5. The results of these calcu- 
lations for a' are plotted in figure 3, which also includes lines of 
constant P/L. Thus, figure 3 can be used to obtain a' and Cq jo deriva- 
tives (Cq/ct being directly related to P/L). Each chart of figure 3 is 

for a constant value of collective pitch, because this parameter is the 
one most likely to remain constant during a change in flight condition 
caused by a gust or maneuver. 


The values of a' and P/L, rather than their derivatives, are 
plotted in figure 3 in order to keep the number of charts down to a 
minimum. Thus, slopes or differences must be measured and some simple 
calculation often must be made in order to obtain the a' and 
C q/o derivatives from the charts of figure 3. 


COMPARISON BETWEEN a' AND aq 


Because of the unavailability of a convenient source of a', most 
stability calculations have heretofore been made on the assumption that 
a' = a^^ because aq was much more readily calculated. The error caused 

by this assumption is shown in figure k, which compares theoretical values 
of a' and aq for p = 0.15, 0.30, and 0.50 and 9 = 0°, 8°, and lii°. 

In general, the comparison between the slopes or increments is quite 
unsatisfactory, and thus it is concluded that the substitution of aq 

for a* frequently yields grossly incorrect rotor stability derivatives. 
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METHOD OF OBTAINING STABILITY DERIVATIVES 


The procedure for obtaining the rates of change of C p/a, a', and 
C Q /o with respect to a, V, Q, and 9 for a helicopter rotor in for- 
ward flight from the charts presented herein is as follows: 

(1) The values of a and R and the trim values of V, p, Cp/cr, 

and P/L must be known. For cases where P/L is not already known, it 
can be obtained from the charts of reference 3. 

(2) The value of collective pitch 9 can be obtained either from 
the charts of reference 3 or by interpolating between the various charts 
of figure 3 presented herein. For a rotor with twisted blades, the col- 
lective pitch on the charts represents approximately the pitch at the 
three-quarter radius. 

(3) The rates of change of thrust-coefficient - solidity ratio with 
rotor angle of attack and collective pitch angle are obtained from fig- 
ure 1(a) where these two derivatives are plotted against tip-speed ratio 
for specified values of rotor solidity. (The upper curves of fig. 1(a) 
are accurate to within U percent and the lower curves are accurate to 
within 7 percent for X from 0.03 to -0.10, Cp from 0.002 to 0.008, 
and a less than 10° in magnitude. For greater accuracy, equations (A3), 
(A10), and (All) of appendix A can be used.) 

(U) The rate of change of thrust-coefficient - solidity ratio with 
tip-speed ratio can be calculated by use of equation (2). The values of 
the five constants in equation (2) can be obtained from figure 1(b) where 
they are plotted against tip-speed ratio for specified values of solidity. 
(The curves of fig. 1(b) are accurate to within U percent for X between 
0.03 and -0.10 and Cp between 0.002 and 0.008. Also equation (2) assumes 

X/p. to be less than 0.3 in magnitude. For greater accuracy, equations (A10) 
and (A13) of appendix A can be used.) 

Derivatives with respect to V and 9. are then obtained as follows: 



and 

8 ( c t/ ct ) = ay. 8 ( C T / a ) 

3ft 9ft 9|JL 


p. 8 ( C t/°) 

9 Op, 


( 10 ) 
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(3) Once the C^/ct derivatives are known, the rates of change of 
longitudinal rotor-vector tilt and torque-coefficient - solidity ratio 
with respect to rotor angle of attack, tip-speed ratio (and, therefore, 
forward or rotational speed), and pitch setting can be obtained for 
various trim combinations of tip-speed ratio, thrust-coefficient - 
solidity ratio, and power input from figure 3, which consists of a 
series of charts for different values of pitch setting. In these charts, 
the longitudinal angle a' between the rotor resultant-force vector and 
the axis of no feathering is plotted against the thrust-coefficient - 
solidity ratio for specified values of tip-speed ratio. Lines of constant- 
power drag-lift ratio P/L, which is directly related to torque coefficient, 
are cross-plotted in the charts. 


(6) The a' and Cq/ct derivatives are obtained from figure 3 by 

measuring slopes or differences and using these quantities in the following 
formulas (sample derivations of these formulas are given in appendix B, 
and a sample problem is worked out in the section following this one): 


da 1 _ 3a 1 _ d( C T/ a ) 

3a da 


/ 3a 1 ) = 1 /3a ' \ 

' const. Cip Jo \ dp / C onst. C. p/a 

/da 1 \ = 1 /da'\ + da' 3 ( C t/°) 

' ^ 'const, a aR W ' const. C T /o 3 ( C t/ c ) 

/ da'~ \ _ p /da 1 \ 

' ^ 'const. Cip/o 'const. C^/cr 


( 12 ) 


(13) 

(lh) 


da 

da) 


const, a 


= _ 1 ± 

a 


da') 

30 / 


at an appropriate combination of p and C^/o. ) 


'da3\ 

3a' 3 ( C t/ ct ) 


(13) 

'const. Crp/o 

3^Ct/ct^ 3p 


Aa ' 

const. C^/ct ^ 



( 16 ) 

are obtained from 

two adjacent charts 

for 
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3a 1 \ = /3a A 


3a' 


3 ( C t/ ct ) 


"/const, a ' "/const. C^/o 3 ( c t/ ct ) 


39 


(17) 


3(C Q /a) 3(C T /a) 

3a * da 


P + 3(P/L) 

L a 


3(C t /ct) 


( 18 ) 


'( c q/ ct ) 


30 


const. Crp/a 


Cm 

= * a 




L 30 J 


const. Cf/cr 


( 19 ) 


3 ( C Q/ a ) 

3(C T /a) 

P ( C T 3(P/L) 

Crp 

3 (P/p) 

L 30 J 

const, a 39 

L a 3(C T /a) 

^ a 

39 J 


>12^1 


dV 


C T /a 

const. Cip/a ^ 


3(P/L) 


P 

L + J 


const. Cip/a 


H c q/°) 

II 

1 

h 

Q 

— - + 11 

3(P/L) 

**\ 

da 

const. Cx /a 

L * 

3m- 

7 

const. C<p/o 


(20) 

( 21 ) 


( 22 ) 


^/r. . /r,\ 

Gijtja 

£ + LL 

3(P/L) 

1 

av J 

const, a ^ 

L ^ 

_ 3m- 

const. C x/a 


>( c T /a) 


ftR 3 m- 


P + 3(P/L) 

L cr 3(Cx/a) 


( 23 ) 
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K g q/°) 


dfl 


J const. 



3(P/L) 





const. C t/ct 


i£ d ( C T/ g ) r P + Ct 3 ( p/ L ) 

ft 3p L ° d(C T /o) 


(24) 


SAMPLE ROTOR LONGITUDINAL-STABILITY-DERIVATIVE CALCULATIONS 


In this section, the longitudinal-stability derivatives, which 
experience has thus far indicated to be significant (except those with 
respect to pitching velocity), are calculated for the rotor for the 
sample helicopter of reference $ by use of the methods indicated in the 
previous section. The procedure for calculating any of the other sta- 
bility derivatives discussed in the previous section are then readily 
deducible and thus are not specifically discussed. For the sample 
helicopter, 


o * 0.07 

R, feet 20 

V, feet per second 80 

p 0.20 

P/L 0.20 

Cr p/CT 0.094 


Determination of P/L and 9 

For cases where P/L is not already available from performance 
calculations or measurements, it can be readily obtained by the method 
illustrated in reference £. If the charts of reference $ are used for 
this purpose, the collective pitch 9 which is needed for some of the 
derivatives can conveniently be determined at the same time. Because 
P/L is already known for this sample problem, 9 can be obtained from 
figure 3. Interpolation between the charts for 9 = 8° and 9 = 10° 
and use of the known values of p, Cf /a, and P/L gives 


9 = 9 . 2 ° = 0.160 radian 

If the sample rotor has twisted blades, this pitch value is the value of 
pitch at approximately the three-quarter radius. 
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Ct /a Derivatives 


From figure 1(a), 


3 ( c t/ ct ) 


3a 


= 0.23 per radian 


3(e T /a) 


30 


O.76 per radian 


From equation (2) and figure 1(b), 

8(C t /o) = _ 3>ii(0<l60 ) + 6.2(0.09U) - 3.1(100)(0.09h)(0.07)(0.l6) 2 

3p 

0.7(1000)(0.09U) 2 (0.07)(0.160) - 3-6(l00)(0.09U) 3 (0.07) 


3( c T /a) 


3p 


= 0.0U 


SIR - — = J4OO feet per second 


SI = = 20 radians per second 

20 


From equations (9) and (10), 

d ( C T/ CT ) = Q • Qk = 0.0001 per foot per second 


3V 


Uoo 


d.fo?/. q .I = (- (O.OU) = -0.000U per radian per second 
3fl V 20 /V ' 


a' and Cq/o Derivatives 


^§1 derivative.- Interpolation between the charts for 9 
3a 

9 = 10° of the slopes of the lines for p = 0.20 gives 

3a' 


= 8° and 


3(Ct/o) 


= 29° = 0.91 radian 
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Using equation (11) yields 

= (0.5l)(0.23) = 0.12 radian per radian 
da 


derivative at constant Qy/o.- Interpolation between the charts 

for 9=8° and 9 = 10° of the increment in a' from p = 0.15 to 

, c T 

[L * 0.25 at — = CL09U gives 


const. C^/ct 


22° = 0.38 radian 


Using equation (12) yields 


f^L-J = 9* ^ = 0.00095 radian per foot per second 

^ const. C T/ a ^ 00 


3a' 

3V 


derivative at constant a.- Using equation (13) yields 


^ 3a' \ 

\ ' const, a 


= Ij^[o.38 + (0.5l) (0.0li)J = 0.0010 radian per foot per second 


9 (C q/o) 

- A v — L derivative.- Interpolation between the charts for 9 = 8° 
9 ^ Q 0 

and 9 = 10° of the increment in P/L from — = 0.08U to — = O.lOU 

a a 

along the lines for p = 0.20 gives 


9 (P/L) 
a(C T /a) 


-ll.8 


Using equation (18) yields 


= (0.20) (0.23)1 0.20 + (0.09U) ( — - 8 )~ 1 = -0.012 per radian 
3a - 
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>(Cq/°) 

39 


derivative at constant Comparison between the charts 


for 9=8° and 9 = 10° at p = 0.20 and ^ = 0.09U gives 


3(P/L) 


39 


_ Q-2 , 5 . 3 _ 0.137 , _ 3.3 per radian 
const . Cip jo ^ 


Using equation (19) yields 




ae 


= 0.20(0.0914) (3-3) = 0.062 per radian 


0(C Q /a) 


39 


-*const. Cj/o 


derivative at constant a.- Using equation (20) yields 


3 ( C Q/ a ) 


39 


= 0.20(0.76)^0.20 + 0.09U(-li.8)] + 0.20(0.0914) (3.3) 


-“const, a 


= - 0.038 + 0.062 

= O.O2I4 per radian 


‘ 3 ( C Q M 


39 


in autorotation.- For purposes of a subsequent 


-'const, a 


discussion. 


~ 3 (°Q/ ct ) 


L 39 


for the sample helicopter in autorotation 
const, a q 


Urn p 

4 ! _ . 11 — n OC\ £. = n noli hill £. = n in+.ornnlst.i nn 

JL 4D I1UVV UU Ud,±UCU . r\JL P — T w 3 ~ ~ — 

O’ -b 

between the charts for 9 = I 4 0 and 9=6° gives 

9 = 3.6° 


5(P/L ).,.. 

3 (C t /ct) 


3 (P/L )1 

L J const. C^j/ct 


= — I4 - I4 


= 3.1 per radian 
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Thus, substitution in equation (20) yields 


_ 39 _ 


const, a 


0.20(0.76)[o + 0.09M-U.U)] + 0.20(0.09U)(3.1) 


= -0.063 + 0.058 
= -0.005 per radian 


ROTOR STABILITY CHARACTERISTICS 


Many significant stability characteristics of the helicopter rotor 
may be deduced from the charts presented herein. A number of these 
characteristics are now discussed. As further illustrations of the use 
of the charts, a number of specific stability problems and current sug- 
gestions as to means for improving helicopter stability are also discussed. 


Rotor-Thrust Variations; Load Factors 


Examination of figure 1(a) shows that the variation of 


K C T/ q ) 

3a 


with p is approximately linear (about 1.2 power). Because (at constant 
rotational speed) rotor thrust is proportional to C-j and forward speed 


is proportional to p, the load-factor increase per degree rotor angle- 
of-attack change also increases in approximately linear fashion with 

velocity. This result is in contrast with the variation with which 

applies for the airplane. Also, because a fixed gust velocity results 
in progressively lower angle-of-attack changes for increasing forward 
speeds, a sharp-edged gust of fixed velocity tends to produce about the 
same normal acceleration for all airspeeds, in contrast with a linear 
increase with airspeed for the airplane. 
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is much less affected 


by p than is 


Figure 1(a) also shows that while 

3(C T /°) 

3a ' g 

increasing p. A change in p from 0.15 to 0.50 increases - 

50 percent. Thus, for a sudden increase in collective pitch, the increment 
in thrust per degree pitch increase at p = 0.50 is about 50 percent 
greater than that for p = 0.l5- 


, an appreciable increase does occur with 

~ ( C T/ q ) 
39 


about 
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Rotor Longitudinal Stability Characteristics; Variations of a' 

For purposes of this discussion, the following assumptions are made: 

(1) The helicopter being considered is a single-rotor machine with 
flapping hinges on the rotor shaft; therefore, its rotor force vector 
passes through its rotor hub. 

(2) The helicopter center of gravity is in its normal position below 
the rotor hub and lies on the trim position of the rotor resultant-force 
vector. 

For these conditions, pitching moments about the helicopter center of 
gravity depend primarily upon variations in a' and an increase in a' 
always causes a nose-up moment. (For those cases where the center of 
gravity is offset from the trim position of the rotor force vector, thrust 
variations produce additional pitching moments. Offset flapping hinges 
also result in additional pitching moments when there is longitudinal 
rotor blade flapping with respect to the rotor shaft.) 

Power-on flight; variations of a 1 with angle of attack .- The charts 
of figure 3 indicate that, during power-on flight (9 * U 0 to II 4 0 ), the 
helicopter rotor is statically unstable with changes in angle of attack 
(or Cj/a) at constant pitch, rotor speed, and forward speed. This insta- 
bility can be deduced from the positive slope of the lines for constant 
tip-speed ratio. When the helicopter noses up, a increases and thus 
Ct/c? increases. Inasmuch as the slope of the lines for constant p is 
positive, an increase in a' results. Because a' is positive when the 
rotor force vector is tilted behind the axis of no feathering, an increase 
in a' (with no change in longitudinal stick position) results in a nose- 
up pitching moment on the helicopter; therefore, the rotor is statically 
unstable with change in angle of attack at constant pitch, rotor speed, 
and forward speed. 

Comparison of the slopes of the various lines for constant p shows 

that — . increases (that is, becomes more unstable) approximately 
3 (C T /ct) ' “ 

linearly with increasing p at constant Crp jc and 0. Thus, inasmuch 

a(Cr/a) -1 p 

as —L varies approximately with p- 1 -* , Ba'/da varies approxi- 

oa 

mately with p^*^ at constant C^/o and 9. Comparison of the various 

charts of figure 3 reveals that the rotor instability with angle of attack 
increases with increasing 9 and, therefore, in general, with increasing 
rotor-shaft power. This increase in instability can be deduced from the 
increasing slope of the lines for constant p with increasing 9. Thus, 
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the variation of da 1 /da with increases in level-flight forward speed 
(beyond the speed for minimum power) has an exponent even higher 
than 2.2 because of the increase in 0 with increase in level-flight 
speed. The increase in thrust that accompanies an increase in angle of 
attack further adds to the rotor static instability. (See reference 11. ) 
The rapid increase in instability at the higher forward speeds as shown 
experimentally in figure 1; of reference 2 is thus seen to be explainable 
by means of these charts. 


Methods for offsetting unstable variation of a' with angle of 
attack .- Inasmuch as the helicopter rotor is statically unstable with 
angle of attack in forward flight, some means must be incorporated on 
the helicopter to offset this characteristic. One means is to use a 
horizontal tail surface as discussed in references 3 and 11. By use of 
the charts presented herein, the unstable rotor pitching moment per degree 
change in angle of attack can be obtained and used to determine the size 
of tail surface needed to stabilize the helicopter. Another method 
which has received attention is to design the hub linkage in such a way as 
to cause the collective pitch to decrease with increasing coning angle. 
Theory indicates that a-^ decreases with decreasing collective pitch at 

constant C ^ja and p.. (See fig. U») As a result, when a-^ is assumed 

equal to a', analysis indicates that a helicopter rotor with sufficient 
"pitch-cone'* change would be stable with angle of attack at constant rota- 
tional speed, forward speed, and collective-pitch lever position. 


By means of the charts of figure 3, this method may be examined with- 
out relying on the assumption that a^ is equal to a' . Comparison of the 


various charts of figure 3 shows that, for typical cruising conditions for 
current helicopters, improvements in stability should be obtainable. 

Further comparison shows, however, that a' , unlike a]_ does not neces- 
sarily decrease with decreasing pitch at constant C^/cr and p. and 

A 

that at high values of the ratio — — , a ! actually increases with 

C<p jo 

decreasing pitch. Under these circumstances, if a given angle-of-attack 
change is considered, then the pitch-cone change will, by minimizing the 
increase in rotor thrust (and, therefore, Crp/cr), reduce the increase in a ? 


caused by the increased rotor angle of attack. The pitch-cone change can- 
not, however, result in stable moments because some increase in thrust 
must be permitted. Also, if a given increase in thrust is required, as 
for a specific maneuver, under these circumstances the result is a greater 
unstable moment than without the pitch-cone change. The value of the 


parameter ■ ■ ■ 
C^/ct 

c ons t ant and 


above which a 1 increases 
C<p/cr, varies from about 3 


with decreasing pitch at 
at (jl = 0.15 to about 2^- 


at [i = 0.50. Because stability problems appear likely to be more acute 
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for high-speed (over 15>0 mph), high-powered helicopters, additional cal- 
culations were made to determine the probable values of for such 

designs. It was found that such helicopters, if of conventional basic 

^ 0 
design, will necessarily operate at values of - — — higher than 3; 

Cj /a 

therefore, the preceding discussion appears applicable to this type of 
helicopter. 


Power-on fli ght; variations of a 1 with speed .- The charts also 
reveal that, during power-on flight, the rotor is normally statically 
stable with variations in speed at constant pitch, rotor speed, and 
thrust. This stability is indicated by the increase in a' with 
increasing p. at constant Cij/ct. Thus, an increase in forward speed 

results in a rearward tilt of the rotor resultant force, tending to oppose 
the increase in speed. This tilt is, by definition, a stable variation; 
therefore, the rotor is statically stable with changes in speed. The 
spacing of the lines for constant p indicates that the stability with 
speed is approximately constant for various speeds at constant Cx/o 
and 9. However, at very low values of Cip/a and high values of p, for 
values of 9 of 6° and higher, the charts indicate that the rotor becomes 
unstable with speed at constant C f/o. This instability indicates a 

problem for the high-speed helicopter, which will necessarily operate at 
low values of Cr ^/o in order to avoid stalling. 


Power-off flight; variations of a' with speed.- For power-off flight, 
the rotor speed varies during a gradual change in flight condition in such 
a way as to maintain zero aerodynamic torque; although for rapid changes 
in flight condition, rotor inertia prevents the rotor from reaching its 
equilibrium speed for a significant amount of time. The derivatives for 
slow changes in flight condition are thus obtained by following the lines 


for 7=0. 

J-j 


For values of 9 ^ 2°, the autorotating rotor is, unlike the 


power-on case, neutrally stable with speed at constant angle of attack. 
(The case of 0° pitch is discussed subsequently . ) This statement can be 

p 

understood when it is realized that at = 0 and constant 9 and a, 


p and, therefore, a 1 must remain constant. (One point on the charts 
represents the whole speed range.) The physical meaning is that a slow 
variation in forward speed at constant angle of attack results in a pro- 
portional variation in rotor speed, such that the tip-speed ratio and, 
therefore, a* remain constant. Thus, the stability with speed at con- 
stant a is neutral. Inasmuch as C T /cr also remains constant during a 


variation in forward speed at constant a, the accompanying variation in 
rotor speed results in a variation in magnitude of the rotor force. If 
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/ a (c t /ct) \ 

a land because ^ — is always positive, C t/ctJ is reduced with 

increase in forward speed in order to maintain constant thrust (as would 
be done in flight in order to maintain lg normal acceleration), the 
rotor speed increases less than in proportion to the forward speed. Thus, 
p and a* increase and the rotor is stable with speed. These state- 
ments can be verified by noting that a reduction in Ct/ct along a line 

p 

for — = 0 results in an increase in a' and an increase in p. Thus, 

L 

the autorotating rotor at constant thrust is stable with speed as for the 
power-on case; whereas at constant a it has neutral stability with 
speed. 


It might be pointed out that the fixed-wing airplane has these 
same two types of stability with speed. If Mach number and power effects 
are ignored, the fixed-wing airplane has neutral stability with speed at 
constant angle of attack and is stable with speed at constant lift. 


Power-off flight; variations of a 1 with angle of attack .- For the 
same pitch range, the autorotating rotor is, in contrast to the power-on 
case, statically stable with angle of attack (or Ct/ct) at constant forward 

speed. This statement can be deduced by again examining the lines for 

p 

= 0. An increase in a results in an increase in Ct/ct and therefore 

P 

a reduction in p and a' because the lines for ^ = 0 slope down to 


the right. The physical meaning is that, if a is slowly increased, 
the rotor speed will increase and p will be reduced at constant for- 
ward speed. Owing to the stability with changes in tip-speed ratio, 
the reduction in p produces a decrease in a' which overcomes the 

1 due to the increase in a at constant p 


increase m 


(near £ - o). 


and thus a' is reduced. Thus, because the rotor speed varies during 
a slow change in a, the autorotating rotor is stable with slow changes 
in angle of attack. 


Power-off flight; 8 = 0° .- For the case of 9 = 0°, the same ideas 

about variations in a' with speed and angle of attack apply as for the 
higher pitch values except that the variation of Ct/ct with p (and a) 
is negligible. 

Rotor-speed effect; power-on flight.- A rotor-speed effect occurs 

in power-on flight similar to that in autorotation in that if a is 
slowly increased, the rotor speed (and engine speed) will increase some- 
what due to the reduction in Cq/ct (at constant p) . The resulting 

reduction in p reduces somewhat the rotor instability with angle of 
attack. 
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Variations of longitudinal stick position for trim.- These ideas 

on the variation of a 1 with speed and angle of attack can be applied 
directly to variations of helicopter longitudinal-control position 
during steady unyawed flight if it is realized that longitudinal cyclic 
pitch is equal to a' when the flapping hinges and the center of gravity 
are on the rotor shaft and when fuselage pitching moments are negligible. 


Rotor-Torque Variations 


The variation of Cq/ct with 0 at constant p, and a is a sig- 
nificant derivative for the coaxial helicopter which uses differential 
torque obtained by differential collective pitch for yawing control. 

The sample calculations of the previous section confirm results, which 
are understood to have been obtained in flight tests, indicating that 
this means of yawing control involves large variations in effectiveness 
with changes in power. For the sample helicopter in autorotation at 




p. * 0.20 and — = 0.09ll 3 the variation of Cq^ct with 9 at constant 


and a is about 20 percent of the level-flight value and is reversed 
in sign. Thus, a coaxial helicopter relying solely on this type of 
control would have reversed yawing control at this flight condition. 
At some intermediate power condition, the helicopter would have no 
yawing control. Thus, coaxial helicopters apparently must use some 
additional source of yawing control in the low power and autorotative 
conditions for satisfactory flying qualities. Synchropter-type heli- 
copters that rely entirely on differential collective pitch for yaw 
control are understood to have also encountered this large variation 
in yawing control effectiveness. 




CONCLUSIONS 


Theoretically derived charts have been presented from which the 
longitudinal-stability derivatives of a helicopter rotor in forward 
flight (except those depending on pitching velocity) can be conveniently 
calculated for tip-speed ratios from 0.15 to 0.50. Examination of 
several basic or current problems by means of these charts indicated the 
following conclusions: 

1. Comparison between theoretical values of rates of change of a' 
(longitudinal angle between the rotor force vector and axis of no 
feathering) and a 1 (longitudinal flapping coefficient) in general 
shows large differences. Thus, the use of the assumption that the rotor 
force vector is perpendicular to the tip-path plane and, therefore, that 
a-^ = a' (often resorted to for convenience) can lead to large errors in 

stability calculations. 
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2. The load-factor increase per degree rotor angle-of-attack change 
(at constant rotational speed) increases in approximately linear fashion 
with the forward velocity V. This result is in contrast with the 

O 

variation with V which applies for the airplane. 

3. In power-on flight, the rotor is unstable with angle of attack 
and is ordinarily statically stable with speed. The instability with 
angle of attack increases with increasing trim value of tip-speed ratio 
or of power setting but is reduced somewhat if the rotor speed varies 
due to the accompanying change in aerodynamic torque. 

It. The autorotating rotor is neutrally stable with speed at con- 
stant angle of attack, is stable with speed at constant thrust, and is 
also stable with slow change in angle of attack at constant forward 
speed. The autorotating rotor is different from the powered rotor 
because of the rotor rotational-speed variation with slow change in 
flight condition. 

5. The charts confirm the reduction or reversal of yawing control 
which is understood to be encountered in low-powered flight and auto- 
rotation in coaxial and synchropter-type helicopters when differential 
torque, as obtained by means of differential collective pitch, is 
employed as the only source of yawing control. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., November 2$, 1950 
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APPENDIX A 


DERIVATIONS OF C T /a DERIVATIVES 


By combining equations (6) and (7) of reference li, the following 
equation is obtained (omitting the blade-twist term): 


C^i 

~o 


crt 


3A 


+ 2 


2p- 




1 + (& 


= tan a + 2 ^ 


(Al) 


Derivative of C > jja with respect to 0.- Differentiation of 

equation (Al) with respect to 9 at constant p and a yields the 
following equation: 


but 


3(0 t /o) 

^3,1 + 2 

+ C T, 

X 0(x/p) 
^3,1 p 39 l 

39 

4-(i) ! 'j 

a 

2u 

! * & 

3/2 


>= t 


3,2 


>(x Al) , 8(x/n) 3 ( c t/ ct ) 

3 (C T /o ) 39 


39 


(A2 ) 


(A3) 


Substituting for in equation (A2) and combining terms gives 


3 ( C x/ g ) 


39 


£ + 

a 


Ot-q 


2pul + 


V- 


X £t 3 (x/p) 

ll a 3 ( Cm ln\ 

i _ 1 : -i-i-L-i 


l + ( - 


> = t 


3,2 


(All) 
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Thus, 



where 


(A5) 


X C T d(\/p) 



The expression for C can be reduced to a more convenient form as follows: 
Equation (?) of reference k can be rewritten as 



Differentiating with respect to C at constant p and a gives 
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Thus , 


a (xAO 

3(C t /ct] 


^ 2 -H ♦ (M 2 


Combining equations (A6) and (A9) gives 



2p2 [ 

_1 + 

c = 

C«p 



\ 

1 + 

Calculations indicate 

that, 

for 


(AlO) 


tively independent of variations in X and Cj and thus depends prima- 
rily upon p. Hence, inasmuch as t^ ^ and ^3 2 depend upon p only, 

/ 3 (Cm /a) 

equation (A5) indicates that is primarily a function of p 

30 


and a. Thus, the values of 


3(C T /a' 


plotted in figure 1(a), which were 


calculated using X = -0.0$ and C-j = 0.005, are less than U percent in 
error for X between 0.05 and -0.10 and C>p between 0.002 and 0.008. 

Derivative of C>p /o with respect to a.- Differentiation of 

equation (Al) with respect to a at constant p and 0 and repeating 
the procedure of the previous section yields the following equation: 


3(°t/°) 

3a 


b.i 11 

cos 2 0 /l , 


(All) 


where C is given by equation (A10). As mentioned previously, C is 
primarily a function of p for p >0.15. For values of a less than 
10° in magnitude cos^a can be assumed equal to -unity with no more than 

^ , 3( Cm I o ) 

a 3-percent error. Thus, for p -i 0.15, ■■■ - is primarily a function 

da 

a(c T /o) 

of p and a. Thus, the values of -t- plotted in figure 1(a), which 
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were calculated using X = -0.05, C^, = 0.005, and a = 0, are less than 
7 percent in error for X between 0.05 and -0.10, Cij. between 0.002 
and 0.008, and a less than 10° in magnitude. 

Derivative of C^/g with respect to p .- The derivative of C^/cr 

with respect to p at constant a and 9 can be set up conveniently 
as follows: 


31 

{ c t /°) 

const, a, 9 

8 (C T /o) 

d( c T /cr) 

f— ) 

- 

dp 

dp 

x A 3a 1 

const. X,9 

const. X,9 


(A12) 


Differentiation of equations (6) and (7) of reference i| with respect 
to p at constant X and 9 results in the following equation (omitting 
the blade-twist term): 


3(°t/°) 


L dy. J 


const. X,9 




(A13) 


' da\ 


= cos'-a 


const. X,9 


-L - 

..2 


2p | 


^ + X" 


2(x 2 + 


J/2 


(Alh) 


Substitution of equations (All), (A13), and (Alii) in equation (A12) and 
rearrangement of terms yields the following equation: 



For normal values of X/p, the last two terms of the braced expression 
can be approximated satisfactorily by means of the binomial theorem as 
follows: 




P 3 



2 


rv 

/v» 


(A16) 
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Substituting this approximate equation in equation (Al5), expres- 
sing X in terms of Cj and 9 by means of equation (6) of refer- 
ence 1;, and rearranging terms results in the following equation for 


df Cf/ a ) 

* — L at constant 

dp 


a and 9 ’• 


where 


* ( c t/ ct ) _ 


dp 


= kq0 + k2 — + k-jCipG^ + k^ 


C t 2 0 


Crn 3 

^5? 


-V (A17) 


k _ a at 3,2 
k l " 2 3ll 


a ^2 dt 3,l 


ILL 


2 t 3 q dp / t 3 qoC' 

LL — + 


2^ 


k 0 = 


1 dt 2il + 


t 3,l a 


h,i ^ 


ap 


/2 + t^oCX L [ t 3>1 oC' 
i^a 2p ) \a 2p 


ko = - 


r _3_ 1 2_ 

M,l) , U t ^ -| oC 


2*1. 


a k * ^ 

\a 2p 


k) , = 


Ut 


3,2 


• ¥ 


k^ - - 



The parameter C is given by equation (A10). 
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Although, as discussed previously, the parameter C varies somewhat 
with 1 and C T , calculations indicate that the values of k plotted in 

figure 1(b) using X = -0.05 and = 0.005 are less them U percent in 
error for values of X between 0.05 and -0.10 and values of C>p between 

0.002 and 0.008. Examination of equation (Al6) indicates it to be suf- 
ficiently accurate for values of X/p of smaller magnitude than 0.5. For 
special cases with values of X/p larger in magnitude than 0.5, or with 
values of C-p or X outside the range for reasonable accuracy in the 
values of k, equations (A10) and (Al5) can be used. 
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APPENDIX B 

SAMPLE DERIVATIONS OF EQUATIONS FOR USE WITH FIGURE 3 


In the section entitled "Method of Obtaining Stability Derivatives, " 
simple formulas are given for use in obtaining stability derivatives from 
figure 3. In this appendix, the most significant of these formulas are 
derived. The derivations of the remaining formulas are then readily 
deducible and consequently are not specifically discussed. 

da'/da derivative.- If the rules for differentiation are followed, 
da' /da. can be written as 


da 1 = da' 8 ( C t/ ct ) 
da d^Cp da 

da'/dV derivative at constant Cp/o.- As before, the rules for 
differentiation give 

/daj\ = dp/da_|\ = 1 /da'\ 

' 8 ^ const. Cpyh 8 ^ 8 ^ const. Cp /a ^v 8 ^/ const. Cp/a 

da'/dV derivative at constant a.- The derivative da'/dp at 
constant a is first obtained as follows: 


/ da _ / da'X + da 1 ^(^t/ q ) 

' 8 ^ ^const. a ^ 8 ^/ const. Cp jo 8 ( 8 T/ a ) dp. 

The second term takes account of the effect of any change in Cp jo 

that occurs when p is changed while keeping a constant. Then, as 
before, 

aal\ , 3a' a ( C l/ . j 

• const. C T /a 3 ( C T/°) * _ 


/da'\ = dl±/da J_\ 

\ 8V /const, a 8V V^ /const, a 
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K c Q/ a ) 


3a 


derivative.- The coefficient Cg/a is given as a function 


of P/L by the following expression from reference 5: 


2 _ P til 


a = r -oa 


Simplifying this equation gives 


_ P 

o - L ~o 

Differentiating equation (Bl) gives 


However, 


Therefore, 


] ( C Q/°) 


W = P d ( C T/°) + S 3(P/L) 

da L 3a a 3a 


3(P/L) = 3 (P/L) d ( C t/ ct ) 
3a ^(^T/ a ) 


d ( C 0./ a ) _ , , d ( C T/ a ) 

~ H* 


3a 


3a 


P + £t 3 ( P/L ) 

L ct 3^C-r/a^ 


(Bl) 


derivative at constant C-jy'a.- Differentiation of equation (Bl) 
with respect to 0 at constant p and Crp jo yields 


R c q/ a ) 


39 


-* const. Cfj. /a 


= H 




const. C^/a 


d(G 0 /a) 

derivative at constant a.- Differentiation of equation (Bl) 


with respect to 0 at constant p and a yields 
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3 ( c q/°) 


_ P 5 ( C T/ a ) + C T 

a(p/L) 

ae 

const. 

l ^ ae a 

a 

ae 


-'const, a 


However, 


a(p/D 


a(P/L) 


a( c T /a) 

ae J 

const, a 

ae _ 

const. C >p/a ^(^T/ CT ) 

L ae J 


const. 


a 


The second terra takes account of the effect of a variation in Crp jo 
when 9 is changed while keeping a constant. Thus, 


d ( C Q/ CT ) 

a(c T /a) 

p ( c t a(p/L) 

Ct 

a(p/L) 

const. Grj>/o 

ae 

. * ae 

const, a 

L a jo) 

* o 

L 39 J 
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Figure 1.- Concluded, 
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Rotor resultant- 



Figure 2.- Schematic side view of helicopter rotor in forward flight. 


Longitudinal angle between rotor force vector and axis of no feathering , < 7 ' deg 
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(a) 9 = 0 °. 

Figure 3 .- Charts for determining a’ and Cq Jo derivatives. 




Longitudinal angle between rotor force vector and axis of no feathering , q', deg 
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Thrust-coefficient — solidity ratio , Cj/c 


(b) e = 2°. 
Figure 3-- Continued. 



Longitudinal angle between rotor force vector and axis of no feathering a', deg 
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Thrust -coefficient — solidity ratio , Ct/cT 
(c) e = k°. 


Figure 3*- Continued. 
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Thrust-coefficient — solidity rat iOj Cj/<5 


(d) e = 6°. 

Figure Continued. 




Long/tudwa/ angle. beween rotor force vector and axis of no feathering 
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Figure 3.- Continued 




Longitudinal angle between rotor force rector and axis of no feat her mg / a', deg 
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Thrust-coefficient — solidity ratio , Cr/cr 

(g) e = 12°. 


Figure 3* - Continued 



Longitudinal ang/e between rotor force vector and axis of no feathering , a', deg 
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Thrust- coefficient — solidity ratio , Cr/tf 

(h) 0 = 14°. 


Figure 3.- Concluded. 






